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Multiplicities in semi-inclusive deep-inelastic scattering are presented for each charge state of tt* 
and mesons. The data were collected by the HERMES experiment at the HERA storage ring 
using 27.6 GeV electron and positron beams incident on a hydrogen or deuterium gas target. The 
results are presented as a function of the kinematic quantities x^, Q^, z, and Ph±- They represent 
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a unique data set for identified hadrons that will significantly enhance our understanding of the 
fragmentation of quarks into final-state hadrons in deep-inelastic scattering. 



I. INTRODUCTION 

Isolated quarks have never been observed in nature. 
When a quark or antiquark is ejected from a bound 
state of quarks and gluons by absorption of a high-energy 
photon, as it separates from the ensemble, a shower or 
"jet" of hadrons is produced. This process may be con- 
sidered to proceed through the generation of additional 
quark-antiquark pairs from the color-field, which even- 
tually combine with the original quark or antiquark and 
with each other until a configuration of color-singlet mul- 
tiquark states is reached. Understanding this hadroniza- 
tion process is an essential element of a complete picture 
of the interaction of quarks in Quantum ChromoDynam- 
ics (QCD), and is basic to an understanding of the dy- 
namics of quark-quark, gluon-gluon, and quark-gluon in- 
teractions. It is described by the polarization-averaged 
fragmentation function (FF) Dp which is the number 
density of hadron type h produced by the fragmentation 
of a struck quark/antiquark of flavor /. 

The flavor dependence of fragmentation functions pro- 
vides a powerful tool for probing the flavor structure 
of the nucleon in hadron-induced hard scattering and 
in semi-inclusive deep-inelastic scattering (SIDIS). In 
the framework of perturbative QCD (pQCD) at leading 
twist, SIDIS is viewed as the hard scattering of a lepton 
off a quark or antiquark, which subsequently hadronizes 
into, e.g., a final-state pion, kaon or proton (or their an- 
tiparticles) . According to factorization theorems [1, 2], 
SIDIS can be described in leading-twist pQCD by three 
components: parton distribution functions (PDFs), hard 
scattering cross sections, and FFs. The hard scatter- 
ing cross section is calculable from pQCD. The PDFs 
parameterize the flavor structure of the initial hadron 
state. Both the PDFs and the FFs are non-perturbative 
quantities, but in the collinear framework where they are 
integrated over parton transverse momentum, they are 
believed to be universal, i.e., not to depend on the par- 
ticular type of process from which they were determined 
[3, 4]. The data available [5] demonstrate that fragmen- 
tation of a quark (antiquark) of a specific flavor is favored 
for a final-state hadron that contains a quark (antiquark) 
of that flavor as a valence quark (antiquark) . The strong 
flavor correlation is reflected in the magnitude of FFs for 
"favored" and "unfavored" fragmentation. This correla- 
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tion is exploited in SIDIS experiments to probe the flavor 
structure of the nucleon through the technique of flavor 
tagging, e.g., to extract the flavor dependence of quark- 
helicity distributions in the proton and deuteron [6]. It 
can be used to extract quantities of interest such as the 
flavor asymmetry in the light quark sea from meson yields 
in SIDIS on the proton and neutron [7, 8]. 

While the knowledge of PDFs is highly developed, the 
data available to date for FFs have been much more lim- 
ited, particularly for unfavored fragmentation. Most ex- 
tractions of FFs [9-12] rely on high precision information 
from electron-positron annihilation into charged hadrons, 
which is available at high energy from experiments at 
SLAC and LEP (e.g., Refs. [13-16]), and has the ad- 
vantage of not being convoluted with PDFs. However, 
these data do not distinguish between quark and anti- 
quark contributions because they are only sensitive to 
the charge sum of specific hadron species (e.g., 7r+ -I-tt^). 
In addition, most data are taken at the mass scale of 
the Z boson, at which electroweak couplings become 
approximately equal and thus only flavor singlet com- 
binations of FFs can be determined. Also, because all 
available data is at the same energy scale, determina- 
tion of the evolution with respect to the four-momentum 
transfer of FFs is difficult. The database has been 
expanded in recent years with the inclusion of results 
for inclusive single-hadron production in proton-proton 
collisions at RHIC. These include measurements of the 
transverse-momentum spectra of neutral pions at central 
rapidities at PHENIX [17], and at forward rapidities with 
STAR [18], as well as similar measurements for identified 
pions and kaons for forward rapidities at BRAHMS [19]. 

Accurate measurement of normalized yields of specific 
hadrons in the final state in SIDIS, i.e., particle mul- 
tiplicities, provides a means of extracting FFs at much 
lower energy scales than those of the collider measure- 
ments. The HERMES experiment with its highly devel- 
oped particle identification and pure gas targets is ideally 
suited for such measurements. The data presented here 
were extracted from measurements of leptoproduction of 
pseudo-scalar mesons in SIDIS that used the 27.6 GeV 
lepton beam of the HERA storage ring at DESY, which 
operated with electrons or positrons. The extraction of 
multiplicities of pions and kaons separately for positive 
and negative charge provides sensitivity to the individ- 
ual quark and antiquark flavors in the fragmentation pro- 
cess. The data presented here for proton and deuteron 
targets are the most precise results for multiplicities cur- 
rently available at this energy scale. With the inclusion 
of the kinematic dependence of the multiplicities on the 
component of hadron momentum transverse to the mo- 
mentum transfer Ph± , these data reach beyond the stan- 
dard collinear factorization, and access the transverse- 
momentum dependence of the fragmentation process. A 
preliminary version of a subset of these data has already 
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been used in a recent FF extraction [20-22]. The data 
extend those of an earher HERMES pubhcation [23] that 
reported results for pion fragmentation on the proton. 

This paper is organized in the following way: Sect. II 
describes the experimental arrangement, while the anal- 
ysis is detailed in Sect. III. The hadron multiplicities 
are presented in Sect. IV, and in Sect. V they are com- 
pared with leading order (LO) calculations based on re- 
cent global fits of FFs. The results are summarized in 
Sec. VI. 



II. THE EXPERIMENT 

HERMES was a fixed target experiment which used 
the lepton beam of the HERA lepton-proton collider at 
DESY. The lepton ring stored electrons or positrons at an 
energy of 27.6 GeV. Typical initial beam currents were 
between 30-50 mA. The HERMES gaseous target [24] 
was internal to the lepton beam line. It consisted of a 40 
cm long open-ended elliptical storage cell aligned coaxi- 
ally to the beam. The storage cell was made of 75 /im 
thick Al and could be operated with a variety of atomic 
gases. This innovative technique allowed the collection 
of data with little to no dilution from other nuclear ma- 
terial. Part of the data was collected using a polarized 
target generated by an atomic beam source [25], which 
could produce a jet of polarized atomic hydrogen or deu- 
terium with an average nuclear polarization around 85% 
and an areal density of 7.6 x 10^^ nucleons/cm^ (hydro- 
gen) or 2.1 X 10^''^ nucleons/cm (deuterium). The nuclear 
polarization of the atoms was reversed at 1-3 minute in- 
tervals. Unpolarized data were obtained using an unpo- 
larized gas feed system, operating with areal densities of 
lO^'^ to 10^^ nucleons/cm^. The high density unpolar- 
ized target was used only for end-of-fill running where 
the beam current was much lower than the average for 
the full fill. The resulting luminosities were of the order 
of 10^^ to lO'^"^ nucleons/(cm^ s). For this analysis, data 
collected with nuclear-polarized and unpolarized hydro- 
gen and deuterium were used. 

The HERMES spectrometer was a forward angle spec- 
trometer with a geometrical acceptance of ±170 mrad 
horizontally and ±|40- 140| mrad vertically. The HERA 
beam lines passed through the spectrometer and were 
shielded in the magnet area by a horizontal iron plate, di- 
viding the spectrometer into two symmetric halves above 
and below the storage ring plane, thereby limiting the 
vertical acceptance close to the beam line. The spec- 
trometer consisted of a front and a rear part separated 
by a 1.5 Tm dipole magnet. Both parts contained track- 
ing devices, while the back part also contained particle 
identification (PID) detectors. A schematic side view of 
the spectrometer can be seen in Fig. 1; a detailed de- 
scription is available in Ref. [26]. 

Track reconstruction was performed using horizontal- 
drift chambers before (EC 1/2) [27] and behind (BC 1/2, 
3/4) [28] the spectrometer magnet. Combining tracks 
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FIG. 1. Schematic view of the HERMES spectrometer. 

from the front and back part of the detector allowed the 
determination of the particle momenta with an intrinsic 
momentum resolution Ap/p between 0.015 and 0.025 in 
the accessible momentum range. 

The tracks found after the spectrometer magnet were 
also used to identify hits in the PID detectors. Lepton- 
hadron separation was provided by combining the in- 
formation of a Ring-Imaging Cerenkov detector (RICH) 
[29], a preshower detector (H2), a transition-radiation 
detector (TRD) and a lead-glass calorimeter [30]. More 
details on the use and the performance of these PID de- 
tectors can be found in Ref. [6]. 

Charged pions, kaons, protons and anti-protons were 
identified using information from the RICH detector. In 
order to achieve good separation over a momentum range 
of 2 to 15 GeV, the RICH made use of two radiators. 
Particles first traversed a wall of aerogel tiles with index 
of refraction n = 1.03 mounted just behind the entrance 
window. The second radiator consisted of C4F10 gas (n — 
1.0014) filling the volume behind the aerogel tiles up to a 
spherical-mirror array with a radius of curvature of 220 
cm. The array deflected the light produced on to a matrix 
of 1934 photo-multiplier tubes per detector half. 

III. DATA ANALYSIS 

The analysis was performed on data recorded in the 
years 2000, 2004 and 2005. Data taken with both po- 
larized and unpolarized targets are included, where the 
regular switching of the spin direction during polarized 
running results in effectively unpolarized data when com- 
bined. The integrated luminosities for the data sets used 
in this analysis are shown in Tab. I. 

The data were selected using the following criteria: 

• There exists a trigger composed of the signals of 
three hodoscopes (HO, HI, H2) and the electromag- 
netic calorimeter in coincidence with the HERA 
lepton-bunch crossing clock. 

• Data quality criteria are met, including good per- 
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Year 


Hydrogen 


Deuterium 




(Pb'') 


(Pb"') 


2000 


132.7 


167.0 


2004 


31.0 


24.3 


2005 


108.2 


137.3 



TABLE I. Integrated luminosities in pb ^ for the two target 
gases and the three data taking- years included in this analysis. 



formance of the particle identification and tracking 
detectors, and of the data acquisition system. 

• The highest momentum particle identified as lepton 
originates from the target cell and is within the 
geometric acceptance of the spectrometer. 

• All other charged tracks originate from the target 
cell and meet the geometric acceptance conditions. 

The spectrometer acceptance was limited by the 
calorimeter dimensions and by the acceptance of the 
spectrometer magnet. This translates into an acceptance 
in terms of the polar scattering angle of 40 mrad < 9 < 
220mrad with respect to the center of the target. 

A. Kinematic requirements 

The kinematic constraints on events are given by the 
spectrometer aperture and the requirement that the 
event kinematics satisfy the conditions > 1 GeV^ and 
> lOGeV^, which select the deep- inelastic scatter- 
ing (DIS) regime and suppress the region of nucleon res- 
onances, respectively. The relevant kinematic quantities 
are defined in Table II. The limits on the fractional en- 
ergy transfer to the target, 0.1 < ?; < 0.85, exclude the 
region where the momentum resolution starts to degrade 
[26] (lower limit) and the low momentum region (upper 
limit) where the trigger efficiencies have not yet reached 
a plateau as a function of momentum. The upper y limit 
also excludes the kinematic region where large radiative 
corrections to inclusive cross sections are required. 

Hadrons are required to have a momentum of 2 GeV < 
Ph < 15 GeV, given by the operation range of the 
RICH detector. Data that are not explicitly binned in 
the fractional hadron energy z also are constrained to 
0.2 < z < 0.8 to exclude the region at high z, which 
receives a sizeable contribution from exclusive processes, 
and to suppress the target fragmentation region at low 
z. 



B. Particle identification 

Leptons and hadrons are identified based on the com- 
bined responses of the TRD, the preshower detector, the 
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4-momenta of incident and 
scattered lepton V 

4-momentum of the target 
nucleon 

4-momentum of the virtual 
photon 7* 

Energy transfer to the target 

Negative squared 
four-momentum transfer 

Squared invariant mass of the 
photon-nucleon system 

Bjorken scaling variable 

Fractional energy of the virtual 
photon 

Azimuthal angle between the 
lepton scattering plane and the 
hadron production plane 

Fractional energy of hadron h 

Component of the hadron 
momentum, ph, transverse to q 



TABLE II. Kinematic variables in semi-inclusive deep- 
inelastic scattering 

calorimeter, and the RICH detector. This response pro- 
vides lepton-hadron separation with an efficiency larger 
than 98% for leptons with contaminations < 1% and 
99% for hadrons with a lepton contamination of <1%. 
The identification of charged pions and kaons using the 
RICH detector is based on a direct ray tracing algorithm 
that deduces the most probable particle types from the 
event-level hit pattern of Cerenkov photons on the pho- 
tomultiplier matrix. This algorithm is described in the 
Appendix of [31]. 

C. Multiplicities 

The multiplicity Af^ of hadrons of the type h produced 
off a target n is defined as the respective hadron yield 
normalized to the DIS yield. It can be expressed in terms 
of the semi- inclusive cross section d^a^ and the inclusive 
DIS cross section d^crois- 

M^{x^,Qi^,z,PhA^ 

_ 1 f2^ d°Ar''(a:B,Q^^,-Pfe_L,0^) J J, 

dxBdQ^ 

_ 1 ^27r d°a''(3;B,Q^^,-P^_L,0^) 1 . 

— d2„ois(-B>Q^) Jo AxjiAQ^AzAPu^A^,^ 

dxBdQi! 

The Born- level multiplicities, i.e., the multiplicities with 
no limitations in geometric acceptance and corrected for 
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radiative effects and detector resolution, are extracted 
from measured multiplicities binned in 3 dimensions: 
(xb, z, Pii±) when the multiplicities as a function of xb 
are desired, and {Q^, z, Ph±) when they are to be given as 
a function of . Due to the strong correlation of xb and 
in the HERMES data simultaneous binning in these 
two variables would not have an impact on the extraction. 
A 0/i binning has been omitted because of limited sta- 
tistical precision. However, a possible impact of the 0^ 
dependence of the unpolarized semi-inclusive cross sec- 
tion and of the acceptance has been accounted for in the 
systematic uncertainties. The 3-dimensional Born-level 
multiplicites are extracted using correction and unfold- 
ing procedures described later in this section, which take 
into account charge-symmetric background, trigger effi- 
ciencies, exclusive vector meson production, kinematic 
and geometric acceptance effects, and smearing due to 
radiative effects. After all corrections and unfolding pro- 
cedures are applied, the final multidimensional binned 
data are available, and can be integrated to yield the 
multiplicities in a one- or two-dimensional binning in the 
desired variables, e.g., the fractional hadron energy z. 



the PID algorithm. Truncating the X row and inverting 
the matrix yields a relation that gives a weight, the cor- 
responding element of -Ptrunc with which each identified 
hadron is counted as pion, kaon or proton: 

^ ~ -Ptrunc ' ^- (3) 

The uncertainties due to RICH unfolding were estimated 
to be less than 0.5% for pions and less than 1.5% for 
kaons. 

3. Trigger efficiencies 

The required trigger combines information from three 
hodoscopes and the electromagnetic calorimeter. The ef- 
ficiencies of the individual detectors are extracted using 
special calibration triggers, yielding an overall efficiency 
depending on the track momentum and the event topol- 
ogy (e.g., events with one or two tracks) that ranged from 
95% to 99%. The events are weighted with the inverted 
efficiency factor. 



1. Charge-symmetric background 

Isolated high-energy leptons from charge-symmetric 
processes like the tt'^ Dalitz decays and photon conver- 
sion into e+e" pairs can produce a signature indistin- 
guishable from that of DIS events. This background is 
most significant at low and is much stronger for inclu- 
sive scattering than for SIDIS. It is taken into account by 
subtracting from the measured DIS or SIDIS event rate 
the number of corresponding events for which the leading 
lepton has a charge opposite to that of the beam parti- 
cles. The correction to the multiplicities is of the order 
of 1% for DIS and 2% for SIDIS at low z. 



2. RICH unfolding 

As described in Sec. IIIB, the hadron identification is 
based on the response of the RICH detector. The effi- 
ciency of the detector is encoded in 3 x 4 matrices binned 
in momentum, charge, and event topology (number of 
tracks in one detector half). They relate the vector of 
true hadron type Th with the vector of identified hadron 
type h: 
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pP 
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pX 



pK 
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pP 

P 

pX 
P 




(2) 



Here, is the probability that a hadron of true type t is 
identified as hadron of type i. The superscript X refers to 
unidentified hadrons. The matrices are extracted from a 
Monte Carlo simulation of the detector response that uses 



4- Exclusive vector-meson contribution 

Exclusive production of vector mesons (p*^, iv, or (j>) 
can be described in the Vector Meson Dominance (VMD) 
model as the fluctuation of the virtual photon into a qq- 
pair before its interaction with the target nucleon. These 
vector mesons subsequently decay into lighter hadrons 
that are then found in the final hadronic state. The cross 
sections for the exclusive production show a depen- 
dence and can be considered as higher-twist effects. They 
do not involve the fragmentation of quarks originating 
from the target nucleon. If fragmentation functions were 
to be extracted from multiplicities that include such an 
exclusive production, they would be process dependent. 
For this reason the data presented in this paper have been 
corrected for hadrons stemming from these processes, but 
the final tabulation includes data with and without this 
correction. 

The fraction of final-state hadrons originating from ex- 
clusive vector-meson decay was evaluated in each kine- 
matic bin using the PYTHIA Monte Carlo generator. 
This PYTHIA version incorporated a VMD model tuned 
to describe exclusive p° production at HERMES [32]. 
Since PYTHIA can only simulate proton or neutron tar- 
gets, the values for deuterium were constructed as the 
combination of the values for these nucleons. The ma- 
jor contribution due to exclusive vector mesons to the 
final state hadron sample arises in the form of pions orig- 
inating from decay. Due to its anisotropic decay-angle 
distribution, pions from decay are concentrated at low 
and high z. For the low statistics high z region near 
z — 1, it is estimated that up to 50% of the charged 
pions originate from while at low z this fraction is 
negligible due to the much larger yield from other chan- 
nels. For kaons, the contribution from decay is less than 
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FIG. 2. Fraction of mesons generated by the decay of ex- 
clusive vector mesons as a function of z, from PYTHIA (see 
text). The widths of the bands indicate the uncertainty in 
the corresponding fractions. The vertical dashed lines are the 
limits in z used in the multiplicity extractions. 



10%. Both the SIDIS and the DIS yields which determine 
the multiplicities are corrected bin by bin for exclusive 
vector meson production. The correction for SIDIS is 
presented in its z projection in Fig. 2. The correspond- 
ing integrated correction for DIS is 0.041 for the proton 
target and 0.044 for the deuteron target. The bands in 
Fig. 2 for pions describe the range in the correction which 
correspond to a la variation in the parameterization of 
the exclusive p cross section. For the (f) meson, the bands 
are half of the total variation in the correction that re- 
sults from varying the parameterization of the exclusive 
(j) cross section over its total range of uncertainty. 



5. Acceptance and radiative effects 



Other effects that influence the extracted multiplicities 



are 



• QED radiative effects, such as vertex corrections 
and the initial and final state radiation of real pho- 
tons by the incoming or outgoing lepton, that alter 
the hard scattering amplitude and can mask the 
true kinematics at the Ij* vertex, 

• the limited geometric and kinematic acceptance of 
the HERMES spectrometer, 

• the detector resolution. 



These effects were evaluated using two Monte Carlo 
simulations that use the same LEPTO/JETSET [33, 34] 
event generator. The first simulation included QED ra- 
diative effects and a full detector simulation. Radiative 
processes and vertex corrections were simulated using the 
RADGEN generator [35]. The leptons and hadrons pro- 
duced were tracked through a GEANT3 [36] model of 
the HERMES spectrometer. Tracks were reconstructed 
using the same algorithm [26] used for real data. 

As this Monte Carlo simulation contains both the 
generated and reconstructed properties of any recorded 
track, bin to bin migration of hadron h can be quanti- 
fied in a migration matrix n^{i,j), where i (= 1 . . . s) 
refers to the bin number based on reconstructed prop- 
erties and J (= . . . s) to the bin number based on the 
generated (Born) properties. The additional Born index 
j = refers to background events generated outside the 
acceptance that migrate into it. The probabilistic infor- 
mation can be summarized in a smearing matrix 



or 



(4) 



The vector n'^amij) is obtained from the second Monte 
Carlo simulation, generated without radiative or instru- 
mental effects. 

As an example, the SIDIS yield of events without ra- 
diative or instrumental effects, i.e., the Born yield, Y^^^^ 
is related to the measured SIDIS yield Y^oas 



^Born(j) 



E 



(j,z)K_(z)i?-n''(*,0)]. (5) 



Here R = <i'ckcd(0/ i^„?c?.(0 where nf,'^^^,{l)) 
is the number of inclusive DIS events in bin I in the MC 
simulation with radiative and instrumental effects and 
^incas(*) is the measured number of DIS events in bin i. 
The matrix [S'^] ~ {j, i) is the inverted SIDIS smearing 
matrix with j > 0. Because the Born multiplicity is 
^Born/^Born' tcrms of mcasurcd quantities it is given 

by 



rn(j) 



E, [S' 



DISl 



(.?■ fc) [yD<JS(A:)i?-nDiS(A:,0)] 



(6) 



If the kinematic binning were fully differential in all 
kinematic variables, in the limit of infinitely small bins 
the smearing matrix defined by Eq. 4 would be inde- 
pendent of the model embodied in the Monte Carlo 
event generator for j > 0. If the smearing correction 
were applied after yields were integrated over some kine- 
matic variable (s) on which both the SIDIS yields and the 
spectrometer acceptance may depend, the accuracy with 
which the smearing matrix corrects for acceptance effects 
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would depend on the success of the model in describing 
this yield dependence. Hence even for those results pre- 
sented here as a function of only one kinematic variable, 
this correction is applied using kinematic binning in three 
dimensions, which is defined in Tabs. Ill to VI. 



D. Systematic uncertainties 

1. RICH unfolding 

The unfolding algorithm for correcting for the ineffi- 
ciencies of the RICH detector is based on probability ma- 
trices as described in Sec. IIIC2. These probability ma- 
trices are obtained from a Monte Carlo simulation of the 
RICH detector, which produces events that are analysed 
with the same particle identification algorithm (described 
in the Appendix of [31]) as that used to analyse exper- 
imental data. The algorithm is sensitive to the number 
and distribution of background hits on the photomulti- 
plier matrix. Different background scenarios result in a 
set of slightly different probability matrices. The varia- 
tion of the final result when using the different P matrices 
has been assigned as systematic uncertainty. The most 
sensitive region is at low momentum (2-4 GeV). Due to 
their lower relative flux, the uncertainty in the fraction of 
identified kaons generated by other misidentified hadrons 
is much larger than that for pions. For the multiplicities 
presented as a function of 0, this contribution to the sys- 
tematic uncertainty is typically less than 0.5% for pions 
and 1.5% for kaons. 



2. Azimuthal asymmetries in the unpolarized cross section 

The unpolarized cross section contains two terms gen- 
erating a cos((/)h) and cos{2(j)h) modulation in the az- 
imuthal angle 4>h- These two terms arise, in part, from 
asymmetries generated by the Cahn [37-40] and the 
Boer-Mulders [41, 42] effects. The moments of these 
modulations have been extracted at HERMES [43]. The 
effect of these modulations is not included in the Monte 
Carlo simulation used for acceptance corrections. The 
simulation assumes an azimuthally uniform PYTHIA6 
production at the Born level. In order to study the in- 
fluence of these moments, a second extraction was car- 
ried out using a simulation in which the events gener- 
ated at the Born level are reweighted by means of a 4- 
dimensional parameterization of the cosine modulations 
extracted from HERMES data. The difference between 
the results for the extractions with the original and the 
reweighted simulations was taken as the systematic un- 
certainty arising from the effects of azimuthal asymme- 
tries in the unpolarized cross section. For most xq bins 
it is less than 2%. 



3. Year dependence 

Both the proton and deuteron data samples are combi- 
nations of data taken over three different years. Compar- 
ing the results from the individual years led to differences 
larger than expected from statistical fluctuations, which 
could not be assigned to specific sources, but must be at- 
tributed to a combination of differences in the data tak- 
ing conditions (e.g., running with polarized/unpolarized 
targets). The differences are mainly concentrated in the 
low momentum region. The largest deviation in the mul- 
tiplicities between any individual data taking period and 
the overall result has been taken as a systematic uncer- 
tainty. With the exception of a few bins, it is never larger 
than 2% of the measured value. 



4- Finite bin width 

The unfolding algorithm which corrects for acceptance, 
limited detector resolution and radiative effects is inde- 
pendent of the Monte Carlo model used to extract the 
smearing matrix if the kinematic bins are fully differ- 
ential and of vanishing width. For this case the effects 
of imperfections in the model on which the SIDIS Monte 
Carlo is based are negligible. The potential impact of any 
residual model dependence due to the finite bin width 
is estimated by scanning over the JETSET parameters 
which control the Monte Carlo event generator [44] . This 
scan is carried out in the space in of the fit to mea- 
sured multiplicities. Nine parameters which constrain 
various features of the fragmentation process are varied 
in the scan. An eigenvector basis approach to the Hessian 
method [45] is used to generate nine-parameter vectors 
which are uncorrelated orthogonal combinations of the 
input parameters to the scan. The intersections of these 
eigenvectors with the contour which lies 68% above 
the best-fit minimum in the scan space, provides input 
parameter sets that characterize the corresponding un- 
certainties in the multiplicities arising from those in the 
JETSET model parameters. The largest deviation of the 
multiplicities from the values extracted with the stan- 
dard version of the Monte Carlo is taken as a systematic 
uncertainty. This uncertainty does not exceed 3-4%. 

IV. HADRON MULTIPLICITIES 

The Born multiplicities, after correction for exclusive- 
vector-meson production, are presented in Figs. 3 and 7. 
In Fig. 3 they are shown for both proton and deuteron 
targets as a function of the energy fraction z in four pan- 
els corresponding to the type of final-state hadron. The 
individual panels compare data for a given hadron type 
taken with a hydrogen (full circles) or deuterium (empty 
squares) target. Error bars on the points for the statis- 
tical uncertainties are too small to be visible. The sys- 
tematic uncertainties are given by the error bands. The 
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FIG. 3. Multiplicities corrected for exclusive vector mesons 
as a function of z from a hydrogen target (full circles) and a 
deuterium target (empty squares). Error bars for the statis- 
tical uncertainties are too small to be visible. The systematic 
uncertainties are given by the error bands. 



FIG. 4. Comparison of measured Born multiplicities for a 
proton target with and without the correction for exclusive 
vector mesons as a function of z. The open (closed) circles 
include (exclude) exclusive vector meson production. The 
statistical error bars are too small to be visible. 



XB 



0.023 - 0.085 - 0.6 



z O.l - 0.15 - 0.2 - 0.25 - 0.3 - 0.4 - 0.5 - 0.6 - 0.7 ■ 
0.8 - 1.1 

Ph± (GeV) 0.0 - 0.1 - 0.3 - 0.45 - 0.6 - 1.2 



TABLE III. 3D binning used for the unfolding correction of 
those multiplicities presented as a function of z (Figs. 3 and 
5). 



z bins are defined in Tab. Ill (z), which together with 
Tabs. IV {Ph±), V (a;B) and VI (Q^) tabulate the binning 
used in the subsequent multi-dimensional representation 
of the multiphcities presented in Fig. 7. To indicate the 
importance of the correction for exclusive vector meson 
decay, the multiplicities for a proton target of pions and 
kaons versus z with the fraction of mesons coming from 
all processes involving exclusive vector meson decay in- 
cluded (open circles) or subtracted (closed circles) are 
compared in Fig. 4. 

The results shown in Fig. 3 for pions (top row) from a 
proton target are in good agreement with HERMES data 
pubhshed earlier [23]. The proton tt"*" multiplicities are 
slightly larger than those of the deuteron, while for 7r~ 
the reverse is true. The ratio of ir^ to tt^ for the proton 
(deuteron) ranges from 1.2 (1.1) in the first z bin to 2.6 
(1.8) in the last z bin. These results can be attributed 
to the dominance of scattering off the u quarks, reflect- 
ing the fact that the fragmentation process u — Tr~^ is 



favored, while the process u — > tt^ is unfavored. With 
rising z, this effect is enhanced. Similarly, the higher 7r~ 
multiplicities for the deuteron are a result of the increased 
fraction of d quarks in the target and of the stronger fa- 
vored fragmentation to the tt~ from the neutron. The 
K"*" multiplicity for the proton is slightly larger than for 
the deuteron, while within errors for the multiplic- 
ities are equal. The ratio K+/K^ rises from 1.5 to 5.7 
(proton target) and from 1.3 to 4.6 (deuteron target) for 
0.2 < z < 0.8 , reflecting the fact that K~ cannot be pro- 
duced through favored fragmentation from the nucleon 
valence quarks. 

Fig. 5 shows the asymmetry 



Ah deuteron 

^d-p — 



proton 



deuteron ' -""proton 



(7) 



between the hadron production off a proton target and 
a deuteron target. Because of the near equality of the 
proton and neutron multiplicities (see Fig. 3), the mag- 
nitude of this asymmetry is small, but it reflects details 
of the quark structure of the targets. The negative values 
for 7r+ and the positive values for tt" reflect the different 
valence quark content of the target nuclei. The mea- 
sured asymmetry in Fig. 5 is more pronounced in the 
high z region for kaons. For positive kaons it is similar to 
that of TT"*", while the corresponding asymmetry for neg- 
ative kaons is near zero, except at higher z, suggesting 
that negative kaons are less sensitive to the valence quark 
content of the target. A LO calculation (see Sect. V) of 
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FIG. 5. The asymmetry A^_^ as defined in Eq. 7 as a 
function of z, for the multiplicities shown in Fig. 3. The values 
given by a LO calculation using fragmentation functions from 
DSS Ref. [22] and parton distributions from CTEQ6L [46] 
are given by the solid curves. 



FIG. 6. The asymmetry j4jj_p as defined in Eq. 7 as a 
function of a;B, for the pion multiplicities shown in Fig. 7. 
The statistical uncertainty is shown by the error bars, while 
the systematic uncertainty is given by the error bands. 



v4^_p shown in Fig. 5 reproduces the measured values 
for positive charge, but strongly overpredicts the asym- 
metries for negative charge. The same trend of negative 
values for tt"'" and the positive values for tt" is evident in 
Fig. 6 where the asymmetry ^5_p is plotted as a function 
of x-Q for four slices in z. There is no strong dependence 
of ^5_p on the slice in z or on xb- 

When the statistical precision permits, binning of ex- 
tracted multiplicities in two or more dimensions can pro- 
vide useful insights into the correlations between kine- 
matic variables and allow the separation of the effects of 
PDFs and FFs. In Fig. 7 the multiplicities are presented 
for positively and negatively charged pions and kaons as 
a function of transverse hadron momentum P/i_l, xb, and 
for the four slices of z between 0.2 and 0.8. The fea- 
tures of the Phi_ distributions result from the combined 
effects of the initial transverse motion of the struck quark 
in SIDIS, of its emission of soft gluon radiation, and of 
the component of transverse momentum which is gen- 
erated by the fragmentation process. The multiplicities 
measured at HERMES as a function of P^i peak at in- 
creasing P/ij^ as z is increased. For negative kaons the 
distribution versus Phi. is wider than for positive kaons 
and may indicate a more extended fragmentation pro- 
cess. In the LUND model [47], which is also used in 
the HERMES Monte Carlo, fragmentation is modeled by 
string fragmentation in which the color field between the 
quarks is the fragmenting entity. That field is modeled 
as a string with a potential energy that increases with 



XB 0.023 - 0.085 - 0.6 

z 0.1 - 0.2 - 0.3 - 0.4 - 0.6 - 0.8 - 1.1 

Phx (GeV) 0.0 - 0.1 - 0.2 - 0.3 - 0.4 - 0.5 - 0.6 - 0.7 - 0.8 - 1.2 



TABLE IV. 3D binning used for the unfolding correction of 
those multiplicities presented as a function of Phi^ and z 
(Fig. 7). 

increasing quark separation. When the separation is suf- 
ficiently large the string breaks. At each string break 
new quark-antiquark pairs are generated. This typically 
leads to more string breaks for unfavored fragmentation, 
which is expected to be more relevant for negative kaons. 

In the context of the multiplicity for tt"*" and 7r~ be- 
ing shown as a function of xb, it should be noted that 
at HERMES kinematics there is a strong correlation be- 
tween a;B and . For the binning given in Tab. V, the 
average value of rises from 1 (GeV)^ in the lowest xb 



0.023 - 0.04 - 0.055 - 0.075 - 0.1 - 0.14 - 0.2 - 0.3 
- 0.4 - 0.6 

z 0.1 - 0.2 - 0.3 - 0.4 - 0.6 - 0.8 - 1.1 

P/ix (GoV) 0.0 - 0.3 - 0.5 - 0.7 - 1.2 



TABLE V. 3D binning used for the unfolding correction of 
those multiplicities presented as a function of sb and z (Figs. 6 
and 7) . 
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FIG. 7. Multiplicities of pions (left panels) and kaons (right panels) for the proton and the deuteron as a function of Ph±, xb, 
and in four z bins. Positive charge is on the left and negative charge is on the right of each panel. Uncertainties are as in 
Fig. 3. 
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TABLE VL 3D binning used for the unfolding correction 
of those multiplicities presented as a function of and z 
(Fig. 7). 

bin to 10 (GeV)^ in the highest bin. A prominent feature 
of the multiphcity data is their kinematic dependence 
on xb, z, and Q^. These data are measured at a much 
lower energy scale than those of most measurements of 
fragmentation yields. Nevertheless, the same qualitative 
trends observed at higher scales are present in the data. 
A strong dependence of the multiplicities on z and a weak 
or vanishing variation with xb and confirm the contin- 
uing applicability of the standard coUinear factorization 
and the universality of distribution and fragmentation 
functions in the SIDIS cross section at the low scales of 
HERMES, (g2) = 2.5 GeV^ and « 10 GeV^. 

A tabulation of the data is presented [48] in four 3- 
dimensional decompositions corresponding to the bin- 
nings given in Tabs. Ill (z), IV (^^j.), V (xb), and 
VI (Q^). Because of the unfolding procedure used to ex- 
tract the Born multiplicities from the measured values, 
the data in the various kinematic bins are statistically 
correlated. In using the data tabulated, the accompa- 
nying covariance matrix must be considered to insure 
that statistical uncertainties are not overestimated. In 
addition, the results of each of the projections of these 
data discussed in the paper are presented. This complete 
data base has been generated including both multiplici- 
ties which have been corrected for exclusive vector meson 
production, as presented in this paper, and uncorrected 
multiplicities. 

An earlier extraction [7] of kaon multiplicities, used 
in an evaluation of the distribution of strange quarks in 
the nucleon, will be superseded by data extracted us- 
ing the new multiplicities presented here. In particu- 
lar, the previous extraction was accomplished using one- 
dimensional unfolding and multiplicities were presented 
for hadron momenta larger than 2 GeV. A reevaluation of 
the strange quark distribution using the newly obtained 
more accurate kaon data will be presented in a separate 
forthcoming paper. 



V. COMPARISON OF MULTIPLICITIES WITH 
LO CALCULATIONS 

To date, analyses of FFs [9-12] have been carried out 
in the framework of coUinear factorization. In this ap- 
proximation, the multiplicity is defined as the integra- 
tion of Eq. 1 over Ph^- In one such model, the LO QCD- 
improved quark-parton model, the hadron multiplicity as 
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FIG. 8. Comparison of the vector-meson-corrected multiplic- 
ities measured for various hadrons with LO calculations us- 
ing CTEQ6L parton distributions and two compilations (see 
text) of fragmentation functions. Also shown are the values 
obtained from the HERMES Lund Monte Carlo. The statis- 
tical error bars on the experimental points are too small to 
be visible. 

a function of z and is given by 

Nnis(Q^) dz T.fe}IlqfixBm<i^B ' 

where the sum is over quarks and antiquarks of flavor 
/, and 6/ is the quark charge in units of the elementary 
charge. The multiplicities in this LO approximation are 
a reasonable starting point for comparing the HERMES 
results with predictions based on fragmentation functions 
resulting from global QCD analyses of all relevant data. 

A comparison of the multiplicities measured by HER- 
MES for SIDIS on the proton with LO predictions is pre- 
sented in Fig. 8. The multiplicities are calculated from 
Eq. 8 (though integrated only over the accepted range in 
xb of 0.023 to 0.600) using values for the FFs taken from 
two widely used analyses, that of de Florian et al. (DSS) 
[22] and that of Kretzer [9], together with parton dis- 
tributions taken from GTEQ6L. For positively charged 
pions and kaons, the results using FFs from the analysis 
of DSS are in reasonable agreement with the HERMES 
results. For negative charges, the discrepancies between 
data and the results based on FFs from DSS are substan- 
tial, particularly for if ~ where the curve predicted lies 
below the observed multiplicity over most of the mea- 
sured range of z. For 7r~ the results from the DSS anal- 
ysis agree with measurement at low z. For both 7r~ and 
K~ , fragmentation is less affected by u-quark dominance. 
Uncertainties in the less abundant production by strange 
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and anti-u quarks may have a larger impact on the pre- 
dictions than for the positively charged hadrons. Al- 
ternatively, next-to-leading-order (NLO) processes may 
be proportionally more important for 7r~ and particu- 
larly , and the discrepances observed here may signal 
the importance of calculating multiplicities at NLO. For 
kaons the DSS results give a better representation of the 
data than the Kretzer curves. This is to be expected, 
since the DSS analysis included a preliminary version of 
the HERMES proton data in its database. The Kretzer 
results are in substantial disagreement with the multi- 
plicities measured for . In Fig. 8 the multiplicities 
obtained from the HERMES Lund Monte Carlo, in which 
the fragmentation parameters have been tuned for HER- 
MES kinematic conditions [20] are also shown. Inclusion 
of the data reported here in future global analyses should 
result in higher precision in the extraction of EEs, partic- 
ularly those describing less abundant fragmentation pro- 
cesses. 



VI. SUMMARY 

HERMES has measured the multiplicity of charge- 
separated pions and kaons as a function of z, PhjL, 
and produced by SID IS off a hydrogen and a deu- 
terium target. This high statistics data set, which re- 
sult from scattering by pure gas targets of protons and 
deuterons, provides unique information on the fragmen- 
tation of quarks into final state hadrons and will con- 
tribute valuable input for the extraction of fragmentation 
fimctions using QCD fits. The comparison of the results 
from the two targets supports the notion of fragmenta- 
tion into hadrons containing the struck quark as valence 
quark being favored, and will allow more reliable extrac- 



tion of unfavored fragmentation functions. The multi- 
plicities measured as a function of Phi^ will provide con- 
straints on models of the motion of quarks in the nucleon 
in the transverse plane of momentum space, as well as on 
the models of the fragmentation process. 
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